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We find the best-fit cosmological parameters for a scenario of inflation witfi only tlie sufficient 
amount of accelerated expansion for the \4>*' potential. While for the simplest scenario of chaotic 
inflation all observable primordial fluctuations cross the Hubble horizon during the slow-roll epoch, 
for the scenario of just-enough inflation the slow-roll conditions are violated at the largest length 
scales. Performing a numerical mode-by-mode integration for the perturbations on the largest scales 
and comparing the predicted anisotropics of the cosmic microwave background to results from the 
WMAP 7-yr data analysis, we flnd the initial conditions in agreement with current cosmological 
data. In contrast to the simplest chaotic model for the quartic potential, the just-enough inflation 
scenario is not ruled out. Although this scenario naturally gives rise to a modiflcation of the flrst 
multipoles, for a quartic potential it cannot explain the lack of power at the largest angular scales. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

Early Universe Inflation is considered the physical pro- 
cess behind the generation of perturbations that produce 
the large-scale structures observed in the Universe today 
as well as the mechanism that solves the main problems 
of the standard cosmological scenario. As a proposal still 
lacking a fundamental theory from which it could emerge, 
there is a variety of scenarios which achieve the funda- 
mental aim of inflation in different manners. 

Current observational results [P start to constrain the 
application of the chaotic inflation scenario 2J for single- 
field inflation models. A variation of chaotic inflation 
contemplates having only a small amount of accelerated 
expansion, just enough to solve the causality and flatness 
problems of the hot big bang model. In this scenario, the 
initial fast-roll dynamics has observational consequences 
at the larges scales, the physics of smaller scales is deter- 
mined by the subsequent era of slow-roll inflation. This 
setup has been considered to explain the lack of power ob- 
served in the lower multipoles of the CMB [3 Ej . Other 
works have analyzed initial conditions for this scenario 
in the context of a preinflationary radiation-dominated 
Universe 1^ and initial conditions for the perturba- 
tions different from a Bunch-Davies vacuum [7 . 

In previous works O [9] , we proposed a modiflcation of 
the chaotic inflation scenario also with a limited amount 
of exponential expansion which is in accordance with re- 
cent observations. Our motivation considered the fact 
that some theories of fundamental physics, including the 
standard model of particle physics, stop being valid above 
a certain energy scale. In those cases, applying chaotic 
initial conditions would imply a complex potential and 
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a decaying inflaton field. By considering an amount of 
inflation of no more than 50 to 60 e-foldings of expo- 
nential expansion this situation can be avoided. As a 
consequence, the Universe starts out in a kinetic-energy- 
dominated stage with an inflaton potential which is many 
orders of magnitude smaller than the Planck scale to the 
fourth, but is real and well deflned. This situation ini- 
tially causes a violation of the slow-roll conditions at the 
moment of horizon crossing when perturbations are eval- 
uated, but the Universe rapidly joins the attractor (slow- 
roll) regime. 

Regardless of this initial violation of slow- roll, we based 
our previous results on the slow-roll expansion. The aim 
of this work is to formalize the implementation of this 
scenario by integrating the mode equations for inflation- 
ary perturbations in a way consistent with its initial con- 
ditions and flnd those values that best characterize this 
scenario in accordance with current data by means of a 
Monte Carlo integration. Here we limit ourselves to the 
potential only, once the the mode integration is car- 
ried out correctly, one can consider the application of this 
scenario to arbitrary potentials. 



II. MODE INTEGRATION 

For the purposes of our work, we need to integrate 
the mode equations for the cosmological perturbations 
for scalar and tensor modes. The fact that this scenario 
of inflation naturally induces a violation of scale invari- 
ance on observable scales, means that we cannot approx- 
imate the power spectra and higher-order observables by 
a power law, as doing so would rely on the validity of the 
slow-roll approximation. 

Our intention is to obtain the initial conditions that 
predict values for the power-spectra and spectral indices 
which are in accordance with current data. We therefore 
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need to integrate the mode equations numerically and 
insert this new information into the CAMB code [lOj to 
supply it with the spectra and then into the COSMOMC 
code [llj . All relevant information on how this is done is 
provided in the next subsections. 



A. Equations of motion : background 



of the system starts in a stage of kinetic energy dom- 
ination, which later naturally leads to one of potential 
energy domination shortly after inflation starts. Within 
this picture, it is valid to consider an initial condition for 
the function ei^i corresponding to a value bigger than 1, 
[l4] , but for a stage of kinetic energy domination, an ini- 
tial condition close to 3 comes out more naturally from 
Eq. fl, as 02 - > V^. 



As a first approximation, in this work we still assume 
homogeneity, isotropy and flatness from the very begin- 
ning of inflation. The equations of motion for the back- 
ground are then 



3M2 



(1) 



3H,j, + V'=0; 



a prime denotes a derivative with respect to the scalar 
field 0, a dot means derivative with respect to cosmic 
time, H is the Hubble rate, whereas V represents the 
inflaton potential. In this case, we adopt the notation 
V = ^ and Mp = mp/VS^ w 2.4 x 10^^ GeV is the 
reduced Planck mass. 

We perform all integrations with respect to the number 
of e-foldings which means for the differential equa- 
tions that 



dH V 



dN M^H 



3H 



(2) 



dN' 



2V; 



Our convention is as follows: < i?' > 0, then 
the inflaton roll downs the potential from right to left 
and H = dN/dt with t cosmic time, therefore dN > as 
dt > 0. 

In what follows, we make use of the horizon flow func- 
tions e„ [12] , developed as a generalization of the slow-roll 
parameters and defined asPj 



Hi. 
H 



1 dem 

Em dN 



m > 0. 



(3) 



Hi refers to the initial Hubble rate (at = 0). In the 
case of single-field inflation models, the function ei rep- 
resents the ratio of the kinetic energy to the total energy 
density of the field, 



B. Equations of motion : perturbations 

In order to calculate the power spectra, we need to 
integrate the equations that describe the evolution for 
inflationary perturbations. We follow the notation and 
conventions of [15 and write the equations for scalar and 
tensor modes in the following manner: 



dr'^ 



+ k 



z dr"^ 
1 d^a 
a dr^ 



uj^ — scalars. 



itfe = tensors. 



(5) 



where u is defined in jl5] as a gauge-invariant combi- 
nation of field and metric perturbations. The comoving 
wave number k is introduced to calculate the power spec- 
tra through a Fourier expansion. The quantity z is de- 
fined as z = a(j)/H, where a is the scale factor and r is 
conformal time defined as dr = dt/a{t). 

As mentioned before, we are integrating all quantities 
with respect to the number of e-foldings A, in conse- 
quence, the perturbation equations are rewritten and in- 
tegrated as 



d^ul 

dm 



+ 



k 



/s(ei, £2, q) 



dN^ 



(1 



, duu 



k 



Uk =0, 



where /s(ei,e2,e3) = 2 - ei -I- |e2 - |eie2 -I- |e2e3 + and 
friti) = 2 — ei, writing \ ^ and \ ^ in Eq. ^ in terms 
of the horizon flow functions 

When perturbations are generated inside the horizon, 
in the region where k/aH — >■ oo, the modes approach 
plane waves of the form |15j 



+ V 



(4) 



In the scenario of inflation being considered, the ini- 
tial value of the potential is limited to be much less than 
the Planck mass such that : Vi ^ M'^ ^ Afp, where M 
represents a scale that bounds the validity of the the- 
ory under consideration. This means that the evolution 



^ This scheme is completely equivalent to other schemes in writing 
the flow equations |13l . 



Uk{r) 



2k 



(7) 



whereas in the opposite limit, when they are outside the 
horizon ^ — > 0, the solution is — z. 

In our case, in order to set the initial conditions for 
equations ^ , we need to consider that the system is not 
in the slow-roll regime and we are actually starting the 
mode integration from the very beginning of inflation. At 
this point, we make another approximation which needs 
to be improved; we set the initial conditions for the per- 
turbations assuming the Bunch-Davies vacuum Eq. 0. 
This clearly needs to be reconsidered as it should not be 
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necessarily the case that from the very beginning of in- 
flation the part of the Universe on which inflation starts 
corresponds to the vacuum state of the field when there 
is kinetic energy domination. Our initial conditions are 
therefore different from those adopted in 3 for a similar 
situation. 

We also need to establish a way to determine the value 
of the first mode to be integrated, the other modes will 
be generated from this value. We follow the discussion 
of Ref . ( [15] ) , section 11 in which a consistent initial con- 



dition keeps the terms (k 



and (k 



Eq. ([5]) positive at the beginning of inflation. Therefore 
the minimum value of k that we choose is given by the 

maximum of the initial values between the terms 



1 d^z 



and 



1 d^a 
a dr^ 



which are determined by the background ini- 
tial conditions given by this scenario. 

The initial conditions we use for the scalar and tensor 
perturbations for each mode k in terms of N are 



scalars 

_ 1 —ikr 



T 



tensors 

1 —ikr 



dN 



Vk^ 
T 
'Uk 



The expressions we use to calculate the power spectra 
are given by [TO] 



(8) 



T 
Uk 



ZT 



with Zs = CLyfel, zx — a. Once the first mode k^m 
to be integrated is determined as indicated above, we 
consider the mode kmax — 10^ fcmm as the biggest value of 
interest for the mode integration. From them we generate 
200 equally spaced modes to perform the integration of 
equations ^ along with the background equations ([2| 
for each of the modes. 



1. Integration Method 

We followed the structure of the code in Ref. TT] with 
the difference that we do not start the integration at the 
pivot point k — aH but before inflation starts and do 
not Taylor expand any of the quantities being integrated 
as we are working with a specific potential. The struc- 
ture of the mode integration as performed in the code 
was done in the following manner: given an initial value 
of £1 above 1, one has to find the place where infiation 
starts and reset the number of e-foldings to to discard 
the expansion that is noninflationary. At this moment 
one determines the value of the first and last modes that 
will be integrated with the perturbation equations as de- 
scribed in the previous paragraph. Then, the background 
equations are integrated until the end of inflation to find 
the values of the potential, Hubble rate, and quantities 
that will be used to convert each of the modes, whose 



wave numbers are given in GeV to units of Mpc~^ as 
needed for the CAMB code. We assume sudden reheat- 
ing to establish this correspondence : 



^Mpc 



- 1 — , kmin 

h 



;^^exp-^^H*(1.5x lO^*"), (9) 

end 



where c is the speed of light in km/s, h is the Hubble con- 
stant which we take equal to 70, k„iin corresponds to the 
minimum scale in the CAMB code in units of Mpc~^, Ti, 
is the temperature of neutrinos today. Tend the temper- 
ature after reheating, AA^ = Nend — A^*, is the number 
of e-foldings between horizon crossing A^* and the end of 
inflation Nend, H*, the value of the Hubble rate at hori- 
zon crossing and the factor 1.5 x 10'^® is used to convert 
from GeV to Mpc"^ 

In the code, the perturbation equations and the back- 
ground equations are integrated beyond horizon crossing 
until decaying and growing modes of scalars and tensors 
have evolved completely and the amplitudes are frozen- 
in. This happens for tensors at fc « j^aH, a bit later 

than for scalars, since the term ^ in Eq. ^ is bigger 

than the term ^gjf at the beginning of inflation. Con- 
sequently we integrate both spectra up to this value. In 
order to optimize the computational time, we only inte- 
grate the perturbation equations until the Universe joins 
the slow-roll regime and the amplitudes of both spectra 
have a difference of less than 10~^ with respect to the 
corresponding slow-roll expressions. Then we only take 
the later as input into the CAMB code. At the same 
time, the system of background equations is integrated 
to find the moment of horizon crossing at fc = aH and 
use those values for the conversion of units for the scales. 



III. MONTE CARLO INTEGRATION 

One of our main purposes in this work, is to find the 
best-fit cosmological parameters for the scenario of infla- 
tion mentioned in the Introduction by considering initial 
conditions for the background and perturbations as de- 
scribed at the end of subsection III Al and IIIBI Since this 
scenario produces primordial spectra that differ from the 
usual power-law parameterization, we wrote a module 
that can be used to provide the CAMB code with the 
information generated by the mode integration and then 
be introduced, into the COSMOMC code. In our case, 
the information needed to do the mode integration are 
the initial values of the first horizon flow function ei^i, 
the initial value of the scalar field and when varied, 
a range of values for the coefficient A of the potential. If 
this parameter is fixed, then we consider a value of 10~^^ 
to get the amplitude of scalar perturbations correct. 

We reused the parameters in the CAMB and COS- 
MOMC codes which are no longer needed to parameter- 
ize the primordial power spectra and characterize other 
inflationary observables since we determine them directly 
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from the mode integration |^ 

We use WMAP7 data [IJ only to perform the MCMC 
analysis. We vary the following parameters : the physical 
densities of baryons, flbh^ , and dark matter, flDMh^, and 
the reionization optical depth t along with the three pa- 
rameters mentioned before (ei^i, X). We do not vary 
the Sunyaev-Zel'dovich template and fix its value to 1. 
From there the code also varies the ratio of the sound 
horizon to the angular diameter distance 9. The equa- 
tion of state of dark energy is not varied and is set to 
w = -1. The model, as implemented in the COSMOMC 
package delivers the posterior distributions of 13 cosmo- 
logical parameters, but at most 7 of them are indepen- 
dent. That is, one extra parameter compared to a A Cold 
Dark Matter scenario of cosmology implemented without 
varying the Sunyaev-Zel'dovich template and without in- 
cluding tensor modes, with the following parameters : 
flbh'^, flDAih"^, 9, r. Us and Ag. We therefore apply the 
Akaike information criterion to our results to measure 
the difference in goodness of fit. This is presented below 
when we discuss our results. 



A. Initial conditions 

We have mentioned that a natural initial condition for 
the scenario here considered corresponds to the first hori- 
zon flow function being close to 3 which determines the 
initial value for the Hubble rate through the equation 




(10) 



The initial value of the scalar field determines the ini- 
tial value of the potential. From the slow-roll approxima- 
tion, one can give an estimate of the interval where the 
initial value of the scalar field must vary in order to give 
the minimum amount of inflation required to solve the 
horizon and flatness problems of the standard cosmolog- 
ical scenario. Here we consider an interval for the field 
given by (j)i/Mp = [20, 25] and let the data decide how to 
adjust the interval, after each run. 

The code takes into account the fact that, when look- 
ing at different regions in parameter space for ei j and 
(j)i/Mp, there might be a numerical issue concerning the 
way in which the initial value of fc* is chosen; sometimes 
it could be, that the ratio k/aH is too close to 1 already 
from the very beginning of inflation. If one is asking the 
code to determine the moment of horizon crossing and 
does the mode integration as long as this value is less 
than 1, the numerical accuracy might not be enough to 
avoid the evaluation being omitted. The reason is that we 

are choosing kmin — ( ^ ^^t ) which is bigger than the 



corresponding term for the scalars as we want to be as- 
sured that both terms (^P - ^3^) and (jt'^ - \ in 

Eq. ([5]) are positive. If ti^i > 1, kmin is chosen at the mo- 
ment when inflation has just started and ei has become 
just one time step smaller than 1 and fc^i„ ~ a^H^. In a 
situation when the initial ei <§; 1, and the initial k is still 
chosen in the way we do for ei > 1, then k„iin/aH ~ \/2 
which is also very close to 1. This problem does not ap- 
pear in the standard way in which the mode integration 
is done, because the value of k is taken within the range 
reached by observations and the initial value of ei is con- 
sidered to be in the slow-roll regime, then they are not 
correlated in the way ours are. 



B. Regions explored for the initial conditions 



From Eq. (10) one can observe that ei < 3. Therefore 



an initial value of ei could be in principle arbitrarily close 
to 3 from below making the initial Hubble rate approach 
the value of Mp as long as 3 is avoided and the numerical 
accuracy is enough to distinguish them. Apart from this, 
we also explore other possibilities in the Monte Carlo in- 
tegration. We consider values of ei^i as low as 2 with 
the consequence that one has to reduce the interval on 
which the initial value of the field can vary so that the 
potential studied in this case can agree with the data. If 
^ 3, the initial Hubble rate is some orders of magni- 
tude below Mp because the potential in this scenario is 
approximately 8 orders of magnitude smaller than Mp. 

We also explored the interval [0.001,2.999] for the ini- 
tial value of £1. We switched off the lensing for this run, 
as it turned out that such a broad prior on is com- 
putationally intense, as is the lensing itself. Our final 
conclusions are not based on this run, but it provides a 
nice consistency check. 



IV. RESULTS 

We explored several regions and cases in parameter 
space, the most representative ones are summarized in 
Table m For each case we had to determine a new covari- 
ance matrix since we could not use the default covari- 
ance matrix supplied with the COSMOMC program as 
it would overrun the settings for our initial conditions |^ 
We took flat priors on ei^i, 4>i/Mp and A. 

The last column of Table |Tj R — 1, corresponds to the 
Gelman and Rubin R statistic to asses convergence of 
the chains P^. These results are for 6 chains of 200,000 
samples each. We use the features included in the COS- 
MOMC code to propose a good covariance matrix but 
stop updating when the convergence factor R-1 is less 



^ Thanks to Wessel Valkenburg for his advice on this point. 



^ Thanks to Jussi Valiviita for explaining this issue. 
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TABLE I. Summary of main results for Monte Carlo integration showing the allowed ranges for the initial values of ei, the 
field and the quartic coupling A. For all cases 20% of rows have been excluded. Values of ei.i and (pi/Mp outside the indicated 
intervals are disfavored by the data. In all cases the tensor modes have been included. 









A 


lensing 


ind. samples 


burn in steps 


R-1 


1 


[2.85,2.999] 


[23.55,25.4] 


[9.28 X 10"", 1.26 X 10""] 


included 


10 949 


261 


0.0017 


2 


[2.9,2.999] 


[23.7,25] 


[9.28 X 10"", 1.25 X 10""] 


included 


20 058 


188 


0.0023 


3 


[2.6,2.999] 


[22.74,25] 


[9.19 X 10"", 1.27 X 10""] 


included 


7 669 


424 


0.0024 


4 


[0.94,2.44] 


[21.8,25.22] 


[9.33 X 10"", 1.22 X 10""] 


not included 


4 307 


828 


0.0012 



than 0.2 to assure that the chains are strictly Markovian. 
The number of independent samples reported is obtained 
by ignoring 20% of the rows |^ We found that the num- 
ber of independent samples is reduced considerably when 
the BB mode of the polarization is taken into account. 
As there are no useful observational constraints on the 
BB modes available so far, we restricted our analysis to 
runs with TT, TE and EE spectra only. 

In Table |ll] the mean values and standard deviations 
for the Markov Chains described in Table |l] are quoted. 
Their corresponding one-dimensional posterior distribu- 
tions are shown in Figs. [l] and [2] only for the last two cases 
of Table m In all four cases the best fit values of the three 
'noninflationary' parameters flijh^ jft^Mh^ and 9 are in 
concordance with results established by other measure- 
ments (supernovae, clusters, etc.). The optical depth r 
is clearly measured. The same can be said on the self- 
coupling A. The initial conditions of inflation, encoded 
in ci^i and (j>i/Mp are less well constrained. We do find 
good fits to the data (see below), but there is also a de- 
generacy between these two parameters, which results in 
the discrepancy of the marginalized and maximized ID 
posteriors in Figs, [l] to [2] see also Fig. [3] below. 

The results for case 1 in Table |l] correspond to the 
lowest possible value of ei^i and lowest and highest values 
of (pi/Mp in which the quartic potential is compatible 
with the data. A lower value of ei^i or wider interval 
for (j)i/Mp would cause all observable modes to cross the 
horizon during the slow-roll regime and thus would result 
in the exclusion of this potential. Similarly, for the cases 
2 and 3 a higher value of (pi/Mp produces the same effect; 
one starts to recover the slow-roll result for this potential. 

Bigger initial values of the field lead to numerical prob- 
lems when lensing is included. This is not the case when 
we turn off the lensing effect and ti^i and ipi/Mp are al- 
lowed to vary from slow-roll to 2.999 and [20,26], respec- 
tively. One can observe in Fig. ^ that the posterior dis- 
tribution of 4>/Mp exhibits a maximum. For ei^i, there is 
a broad region of almost equally-likely values. This also 
repeats in cases 1 and 3, although for a narrower interval 
and different initial conditions. 

There are some comments concerning the shape of the 
distributions and the number of independent samples 



which are interesting. For £ [2.85,2.999], 4'i/Mp € 
[23.55,25.4], not shown in the table, we obtain very dis- 
torted distributions and observe that this feature repeats 
whenever we explore regions of initial conditions for ei 
below 2.9 and do not include tensor modes. Particularly, 
for the interval ei^i € [2.6, 2.999] without tensors and A 
fixed, also not shown, the posteriors are observed to be 
very distorted, but this is somewhat alleviated by intro- 
ducing the tensor modes and varying A, as can be seen 
from Fig. [T] Reducing the initial value of ei, below the 
favored values in this scenario, causes problems in the 
simulations. 

Regarding the number of independent samples; as 
mentioned, we have adopted for all cases in Table |T] to 
discard 20% of the samples. However, in those cases 
which correspond to ei ^ being almost 3, the burn in of 
the chains occurs earlier. For case 2, discarding only 3% 
of the rows gives 23 813 independent samples, whereas 
considering ei,^ S [2.9,2.999] , <j)^/Mp e [23.7,25], A e 
[9.28 X 10"", 1.25 X 10"i2] with tensors and including 
the BB mode we find only 11 362 independent samples; 
ei,, e [2.9,2.999] , € [23.66,25] and A = 1 x lO"!^ 

including only the TT mode without tensors discarding 
3% of the rows gives 24 328 independent samples. For 
the other distributions, we did not obtain a significant 
improvement in the final number of independent samples 
by repeating the simulation and adjusting of the start 
width and standard deviation estimates in the initial dis- 
tributions, the numbers shown here were somehow the 
best obtained after many attempts. Whether there is a 
physical argument behind this result, that is, whether the 
number of independent samples indicates how physically 
more likely is a prior distribution, or if it indicates the 
region where the scenario is more favored by the data, 
we do not know. 

In Figure|3]the two-dimensional posterior distributions 
of ei.i and (fii/Mp obtained from the Monte Carlo inte- 
gration for case 4 in Table [l] is shown. One can observe 
the degree of degeneracy between these two parameters 
clearly. 



A. Constrained case 



* Thanks to David Parkinson for suggesting this value 



From the results for the distributions of cosmological 
parameters, it is possible to observe that ti i and (pi are 
degenerate with each other. In order to obtain a bet- 
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TABLE II. Summary of main results for Monte Carlo integration. Mean values and standard deviations for examples in Tablejl] 
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DO. DO 




a 


0.00011 


0.005 


0.0021 


0.013 


0.22 


0.022 


0.031 


0.026 


0.067 


0.026 


1.136 


0.0031 


1.86 


3 




0.022 


0.11 


1.037 


0.08 


24.31 


2.77 


-4.53 


0.71 


13.89 


0.29 


10.13 


0.13 


68.75 




a 


0.00032 


0.0054 


0.0027 


0.013 


0.36 


0.102 


0.031 


0.027 


0.08 


0.027 


1.14 


0.0029 


1.92 


4 




0.022 


0.11 


1.037 


0.08 


23.43 


1.69 


-4.53 


0.71 


13.9 


0.29 


10.20 


0.13 


68.83 




a 


0.00032 


0.0054 


0.0022 


0.013 


0.39 


0.43 


0.032 


0.027 


0.08 


0.027 


1.15 


0.0028 


1.94 
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13.7 13.S 13.9 14 14.1 




0.125 0.13 0.135 0.14 



FIG. 1. Posterior distributions of cosmological parameters from case 3 in Table |T] The dotted lines represent mean likelihood 
of samples and solid lines marginalized probabilities. 



ter constraints we have varied another instance: keeping 
(j)i/Mp fixed to a certain value and varying e^^i and A, 
including tensors and lensing and using the TT, TE and 
EE angular spectra, which is represented by the value of 
the parameter {C£s = 3). These results were calculated 
with two chains of 100 000 samples of length each, the 



best-fit values for this case are shown in Fig. |4] and the 
results for a ACDM model without tensors and without 
running is shown in Fig. [Sj 

The posterior distribution produced a flat distribution 
of i and a central value for A = 1.11 x 10~^^±0.03 with 
a convergence R — 1 = 0.0008, an approximate number 
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FIG. 2. Posterior distributions of cosmological parameters from case 4 in Table |l]. 



of 3921 independent samples, and bum in of 40 rows. 

We compare our results with a ACDM model for 6 
parameters: flbh'^, ^oMh"^, 0, t, Ug, Ag giving a total of 
11 parameters along with ujjy, a;„j, z^e, Hq and age of the 
Universe. That is, we have the same number of primary 
parameters in the constrained model and one parameter 
more in the models reported in Table [ij Therefore we 
apply the Akaike information criterion : 



AIC = -21nZ:„,ax + 2A^pa 



(11) 



in order to compare the goodness of fit of ACDM model 
with the four models reported in Table |l] In Eq. 11 C„ 



is the maximum likelihood of the model and iVpar is the 
numb er o f parameters of that model. This is done in 
Table EilB 

The main difference between the ACDM model and our 
results comes from the num ber o f parameters (AA^par — 
1). The last column in Table III shows the value of a 



TABLE III. Akaike information criterion applied to results of 
Table [I with respect to a ACDM model of 6 parameters. 



Model 


In ZImax 


Apar 


AIC 


AAIC 




ACDM 


3737.2 


6 


7486.5 





7474.4 


1 


3737.4 


7 


7488.8 


2.3 


7474.8 


2 


3737.4 


7 


7488.8 


2.4 


7474.8 


3 


3737.4 


7 


7488.8 


2.3 


7474.8 



Thanks to Zong-Kuan Guo for his explanation. 



goodness of fit for these models. Compared to the ACDM 
model, it differs by 0.4 which assures that our scenario 
can fit the data almost as well as ACDM. In terms of 
the AIC, the ACDM with power-law spectrum without 
tensors is slightly favored, but not at a significant level. 

As observed from Figs. |4j [Sj fixing the value of the 
field has the consequence of changing the shape of the 
distributions for the reionization redshift Zre and the op- 
tical depth T. The initial value of ei however seems to 
be allowed to include values lower than 2.9 and we have 
cut the distribution in accordance to the prior imposed. 



FIG. 5. Posterior distributions of cosmological parameters for 
the ACDM model with 6 independent parameters. 



22 22.5 23 23.5 24 24.5 25 



FIG. 3. (color online). Degeneracy between ei,i and (jii/Mp 
for case 4 in Table H] The contours show the 68% and 95% 
confidence limits of the marginalized posterior distributions. 
The shading indicates the mean likelihood of samples. 
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FIG. 4. Posterior distributions of cosmological parameters 
from the prior ei,i G [2.86,2.98], <;/>,/Mp = 25, A G [1.0 x 
10-^^1.23 X 10"^^], Cls = 3 with tensors. 



although the best-fit value is compatible with the argu- 
ment that a natural ei i must be close to 3. 



different scales using their slow-roll expressions expanded 
at first order. We use the following expression to trans- 
late the value of the scale fc* which represents horizon 
crossing from GeV to Mpc~^. This expression was also 
used in [9j for the same purpose and is derived assuming 
sudden reheating. 



kf — 500 exp 



i?*;^(0.002Mpc"'). (12) 



The assumption of sudden reheating is implicit in the 
value of Treh calculated from the temperature at the end 
of inflation. The predictions for the inflationary observ- 
ables and the values of the lower multipoles of the CMB 
for each of the cases presented in Table |l] are presented 
in Table |IV| and |V] respectively. 

In Table |IV| r, is the tensor to scalar ratio, Ug is the 
spectral index of scalar perturbations and dn^/dlnfe the 
running of the spectral index, Nt represents the total 
amount of inflation produced by the model. We present 
the predictions for the observables on different scales: the 
first value corresponds to that on which the Monte Carlo 
integration was done, 0.05 Mpc~^. The second to the 
value that we used in our previous work, 0.01 Mpc~^, 
the third one is that of 0.002 Mpc^^, where the analysis 
of WMAP7 yr is done. N^, are different points on the 
integration trajectories with respect to the number of e- 
foldings where one finds the corresponding value of fc* 
by using Eq. 12 In a previous work |9j, we compared 
the slow-roll trajectories obtained by the WMAP team 
to those of this scenario in the Ug-r plane. We showed 
that the modified initial conditions lead to an increase 
of the spectral index at the largest scales, in such a way 
that allows the A^'* potential to be allowed by current 
data. 



B. Slow-roll predictions 

In this section, we apply the results for the initial val- 
ues of ei^i and ipi/Mp shown in the last section to cal- 
culate the spectra of anisotropics for the CMB and give 
the values for the inflationary observables calculated at 



C. CMB Anisotropies 

We intend to see whether the suppression of power on 
large scales on the primordial power spectra produced by 
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TABLE IV. Values of inflationary observables evaluated at diiTerent scales from the result of the Monte Carlo integration. The 
cases are arranged according to the total amount of inflation produced by each of them. 





Best-fit values: ei (l)i/]\4p^ A 


Nt 




kt {Mpc~'^), TV* 


fc« (Mpc"^), Nt 


fe*(Mpc"^), N, 


3 


2.77, 24.31, 1.1 X 10~^^ 


64.31 




0.05, 59.4 


1.03 X 10"^, 61.0 


2.12 X 10"^, 62.61 








r 


0.27 


0.26 


0.28 








Us 


0.95 


0.95 


1.26 








dus/dlnk 


-9.0 X lO"* 


-8.42 X 10"^ 


-0.9 


4 


1.69, 23.43, 1.1 X 10"^^ 


64.25 




0.055, 59.3 


1.03 X 10"^ 61.0 


2.15 X 10'^ 62.6 








r 


0.27 


0.26 


0.29 








Us 


0.95 


0.95 


1.3 








dUs/dlnk 


-9.0 X lO"* 


-9.77 X 10"^ 


-0.97 



this scenario as seen in |9j could explain the lack of cor- 
relation on large angular scales in the lower multipoles 
of the CMB, [19], [20], [21]. There are other works that 
explore the effects of different inflationary setups to ex- 
plain the lack of power at low multipoles of the CMB: 
using a fitting function and a cut off in a fast-roll era 
[22], a primordial magnetic field [23], anisotropic infla- 
tion [24', using a curvaton model 'SF and isocurvature 
perturbations, 26\. In our work, we do not fit the primor- 
dial power spectra, but allow the numeric solution of the 
mode and background equations to fully determine its 
shape. We use the initial conditions for the background 
dynamics, determined through the Monte Carlo integra- 
tion, to calculate the value of the pivot scale on which 
perturbations are evaluated using their slow-roll expres- 
sions. The assumption of sudden reheating to obtain the 
value of the pivot scale in Mpc~^ gives an uncertainty on 
these numbers. 

There is an important point to mention before pre- 
senting the results of this subsection: We obtained the 
same values for Ctt for all pivot scales fc* considered in 



Table IV This is not the case for the standard power- 
law parameterization of the primordial spectra. In our 
case, the initial conditions for the mode integration are 
not dependent on the value of the pivot scale. We do 
not use this value because the scenario we are investigat- 
ing is not in the slow-roll regime at horizon crossing and 
therefore, setting initial conditions at horizon crossing, 
even if they violate slow-roll do not assure us that we are 
going to recover those initial values which characterize 
this scenario. 

The values of the multipole moments Ci do change 
slightly for different initial conditions of ei and (f>/Mp 
in Table [T] One can observe however, that those values 
are not much lower as compared to the standard power- 
law parameterization and initial inflationary conditions 
in slow roll. In order to be sure that we were giving the 
values of the mode integration correctly into the CAME 
code, we used the power-law parameterization of the pri- 
mordial power spectra in our code and the usual slow-roll 
power-law result of CAMB without any modiflcations. 
For the same initial value of the parameters, we checked 
that the resulting CMB anisotropics for both instances 



at small scales coincided, as was expected. We had to be 
careful not to produce a cut-off in the interpolation since 
it would affect the integration of modes which are deeper 
subhorizon. 

There are some comments to be made in the following 
about how the choice of initial conditions and the way in 
which the primordial power spectra are included into the 
CAMB code. 

In order to include the new primordial spectra to cal- 
culate the Ce values, one must interpolate them so that 
the code can find a wider range of values to the whole 
range of scales reached by observations. We used the 
subroutine SPLINT from numerical recipes to interpo- 
late. It turns out that, for some initial conditions in our 
scenario, this subroutine produces power spectra that are 
negative. We set this output to as long as this happens. 
For those initial conditions, one has at large scales that 
the value of Ctt becomes smaller. Then one can either 
try to fit the quadrupole or octopole separately by trying 
different initial conditions as the initial k changes, or a 
bigger value for the first initial k. Once this is fixed in 
the way explained in Sec. |IIB one can choose to use a 
bigger value for this first scale as it assures that the terms 
multiplied by Uk in Eq. ^ are still positive. We stress 
that this is not the value for the pivot scale at horizon 
crossing used for the Monte Carlo integration but the 
value of the first mode to be integrated at the beginning 
of inflation. 

If one is exploring parameter space for eu and (pi/Mp, 
one has to be careful not to hit regions in which the 
amount of inflation produced is too small. This can either 
happen for a flxed ipi/Mp that gives enough inflation for 
a certain value of ei^i and as the later goes closer and 
closer to 3 the amount of total inflation is reduced, or for 
a fixed ei^i and a small value of the field. Both situations 
cause the interpolation to produce again a negative power 
spectra. 

Having these situations in mind, one can find either ini- 
tial conditions that fit the value of the quadrupole and 
the octopole separately, or increase the value of the initial 
scale to be considered for the mode integration so that 
the final output gives less power in the lower multipoles, 
this would be a complementary way to compensate for 
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TABLE V. Predicted CMB band power ATi = l{£+l)Ce/2TT 
of quadrupole and octopole for all best-fit models shown in 
Table nil 









1 


1 228 


1 139 


2 


1 234 


1 146 


3 


1 228 


1 140 


4 


1 228 


1 140 



.» 




FIG. 6. Initial conditions in concordance with the observed 
CMB quadrupole. 



the situation mentioned at the beginning of this subsec- 
tion. The value of the pivot scale considered in CAMB 
does not have any influence in the mode integration, but 
one can consider varying the value of the first scale to be 
integrated as long as the equations for the perturbations 
are well defined. However, we have not explored further 
this last situation and its consequences for the infiation- 
ary observables and the Monte Carlo integration. 

If one tries the first method, the result is that the up- 
per multipoles are misplaced, this is different from the 
result found in [5 where they use a distribution to fit 
the CMB and different initial values for the scale factor 
considering a preinfiationary radiation dominated epoch. 
Those initial conditions that enter in the first possibility 
are listed in Figure |6] 

By keeping the value of ei.i fixed and finding the value 
of (pi/Mp that gives the desired result for f = 2, we found 
that as long as ti^i > 1, this happened always at a value 
of (j)i/Mp that gave the same amount of total inflation. 
Once ei^i < 1 we did not flnd the same effect. 

By increasing instead the value of the initial scale to 
be considered for the mode integration by a factor of 
60, one can reduce the values of the quadrupole and oc- 
topole arriving to values as low as 346.7 for ^ = 2 and 
344.7 for £ = 3. The values are the lowest ones obtained 
even if one uses bigger scales. In some of those cases the 



interpolation again gives negative values for the output 
of the primordial power spectra and for other values of 
k-minj which could be bigger or smaller, gives a positive 
output throughout the whole interpolation process. The 
quadrupole and octopole could no be fit at the same time. 

For all the models that fit the quadrupole and octopole 
by trying different initial conditions, there is a range of 
N^, or equivalently, of pivot scales in which the infiation- 
ary observables r,ns,dns/dlnk are inside the Icr, 2cr, Ict 
intervals respectively, for WMAP7 with running and ten- 
sors, but not the amplitude of scalar perturbations, which 
turns out to be very small : 2.11 x 10^^ in the best of 
cases, even for those initial conditions with eij < 1. For 
the octopole on the contrary, the amplitude is too big : 
3.31 X 10-^. 



D. Mode Integration results 

Some comments on the results of the mode integration 
are made in this subsection. Wc have mentioned before 
that for a fixed initial value of the field, the bigger ei_i is 
or closer to 3, the less amount of total inflation produced 
by the model there is. If one keeps fixed instead the 
value of the scalar field, in order to have more inflation 
produced, the value of ci^i has to be decreased. That is, 
ei^i and (f>i/Mp are correlated and the only free parameter 
is the total amount of inflation produced Nt- For a flxed 
amount of total inflation, the bigger ei^i is, the bigger the 
initial value of the field must be in order to compensate. 

For the mode integration of perturbations, we fixed 
the initial value of the scale factor to 1 at the begin- 
ning of infiation. As in this scenario there is an epoch 
of noninflationary expansion, we discarded the number 
of e-foldings during that time and reset it to at the 
beginning of infiation. 

In Figs. [7] and [8] we present examples for the primordial 
scalar and tensor power spectra as obtained from a mode- 
by-mode integration and compared to the corresponding 
slow-roll predictions. 

We can observe oscillations in the amplitude of the 
scalar and tensor perturbations. We have seen before 
that the initial value of the first mode to be integrated 
is already very close to horizon crossing due to the way 
it is chosen for this scenario. This means, that we are 
still able to observe the sub horizon oscillations of the 
solutions in Eq. ([s]) before the growing mode dominates. 
In [5], it is said that those oscillations disappear by con- 
sidering an initial condition for the scale factor deeper 
into a previous epoch of radiation domination. We are 
not considering this situation here. However, we have 
checked that if one chooses a bigger initial scale to do 
the mode integration, the cutoff and oscillatory behavior 
observed before disappear completely and the mode inte- 
gration and slow-roll spectra are the same for all scales. 
That is, one obtains pure power-law spectra. 
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FIG. 7. Scalar power spectrum for the initial condition ei,i = 
2.93, ct>i/Mp = 24.67, A = 1 x 10"^^. The main plot shows 
the scales in the x-axis in GeV, the subplot in Mpc~^. 
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V. SUMMARY AND DISCUSSION 

The idea of considering a limited amount of inflation- 
ary expansion as a modification of the chaotic inflation 
scenario implies a choice of initial conditions that have 
observable consequences. In this scenario, the kinetic 
energy dominates before the onset of inflation and conse- 
quently violates the slow-roll conditions at horizon cross- 
ing of observable modes before the system joins the in- 
flationary attractor. 

The best-flt power spectrum of scalar perturbations 
shows a sharp cutoff on large scales. This however, does 
lead to a significant suppression of the angular CMB 
power at the largest angular scales. The reason is that the 
cutoff in a A(/>** model with just enough e-foldings turns 
out to be at scales that are too large to account for a 
lack of power in the CMB at the largest angles. However, 
one can obtain values very close to those of the reported 



quadrupolc and octopolc by considering bigger values of 
the initial scale to be integrated, and also by tuning the 
initial conditions to flnd them separately. This, however, 
causes higher multipoles to have much higher values than 
observed. However these results are only applicable for 
the A(/)^ potential. 

The results of the mode integration for the tensor 
modes show that they have less power, as expected, than 
the scalars and their power spectrum also shows a sharp 
cutoff at large scales. 

An interesting point concerning the initial conditions 
of the inflaton and ei is that the Monte Carlo integration 
gives predictions of best-fit values only for those inter- 
vals specified in the priors. It does not give information 
about the physical significance of the region explored, 
unless that information is contained in the amount of in- 
dependent samples given by the simulation for different 
distributions. We have seen that values ei,, w 3 give 
bigger numbers of independent samples. 

The Monte Carlo integration confirms in the sense 
mentioned above, that a natural initial condition for this 
scenario corresponds to the first horizon flow function be- 
ing close to 3. We have already seen that an initial epoch 
of kinetic energy domination naturally leads to this initial 
condition. 

We have been able to prove with the results obtained 
by the Monte Carlo integration, that the A(/)^ potential 
is inside current observational constraints for single-field 
models of inflation. 

The implementation of this scenario still lacks the con- 
sideration that the initial conditions for the perturba- 
tions cannot be in the Bunch-Davies vacuum and that 
the equations of motion for the background are approxi- 
mated by a flat homogeneous Universe. Both should be 
reconsidered in order to have consistency with the fact 
that there is kinetic energy domination at the beginning 
of inflation. However, as a first approach, we could ob- 
tain enough information of the main behavior for the in- 
fiationary observables and the predictions for the CMB 
anisotropics given by this scenario. One should be able 
now to apply it to different potentials and explore more 
general situations and their consequences. 
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